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PREDICTION OF FRICTION AND HEAT -TRANSFER COEFFICIENTS 
WITH LARGE VARIATIONS IN FLUID PROPERTIES 
by Maynard F. Taylor 
Lewis Research Center 


SUMMARY 


The conventional methods of predicting single- phase turbulent heat-transfer and 
friction coefficients often give values that are in poor agreement with measured values 
when large variations in the fluid properties are present. Heat-transfer and friction 
coefficients calculated by recently reported prediction equations were compared with 
measured values for hydrogen, helium, nitrogen, air, and carbon dioxide from several 
different investigations. 

For the case of heat extraction from a gas the heat-transfer coefficient can be accu 
rately predicted with the relation 

Nu b = 0.023 Re b >8 Pr b ' 4 


and for the case of heat addition with the relation 


Nu b = 0.023 Re b * 8 Pr b ‘ 4 ' T 



where Nu b is the bulk Nusselt number, Re b is the bulk Reynolds number, Pr b is the 
bulk Prandtl number, T is the surface temperature, and T h is the bulk temperature. 

S 

For both heat extraction and heat addition to a gas the friction coefficient can be 
accurately predicted by the relation 


f 

2 



where f is the friction coefficient, Re o is the modified surface Reynolds number, T 

S 

is the surface temperature, and T b is the bulk temperature. 

For both friction and heat transfer the exponent of the temperature ratio T g /T b is 

C 9 = 0. 57 - 
Z x/D 


where x is the distance from the entrance of the test section and D is the inside di- 
ameter of the test section. These prediction equations are applicable to gases flowing 
subsonieally through a smooth tube. 



INTRODUCTION 


Large variations in the physical properties of gases flowing turbulently through 
tubes have been found to greatly affect the heat-transfer and friction coefficients. The 
conventional "reference temperature" methods (i. e. , evaluating the properties and den- 
sity at a reference temperature T x = X(T g - T fe ) + T fe ) can predict coefficients that are 
greatly in error. For example, figure 1 shows that the measured friction coefficient 
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Modified film Reynolds number Re f = VD/v^ 

Figure 1. - Variation of local and average friction coefficients with modified film Reynolds number. Viscosity and density evaluated at film temperature. 

can be as much as three times the value predicted by the Karman-Nikuradse equation 
with the physical properties evaluated at the film temperature T f = (T g + T b )/2: 



= 2 log Re^ 



0 . 8 


(1) 


Errors of the same magnitude are found when measured heat-transfer coefficients are 
compared with those calculated by 
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( 2 ) 


Nu f = 0.021 Re®* 8 Pr®* 4 

Equations (1) and (2) are applicable only for developed flow. 

The use of the ratio of surface to fluid-bulk temperatures T /T, raised to a power 

has been very successful in correlating both friction coefficients (ref. 1) and heat- 

transfer coefficients (ref. 2) with large variations in the fluid properties. The study of 

available friction data (ref. 1) resulted in the recommendation that the modified surface 

Reynolds number Re = VD / v be used in correlating the friction coefficients. It ap- 

s s 

peared that a modified surface Reynolds number of 3000 was critical, and the following 
correlation equations were recommended: 

For Re < 3000: 
s 


2 Re g 

For Re > 3000: 
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= 10.0007 


\ /rrt \-0. 5 

0. 062 5\ ( s \ 

Re®* 32 )\ T b/ 
s / x ' 
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( 4 ) 


The recommended correlation equation for heat addition (Tg/T^ > 1. 0) to hydrogen, 
helium, and nitrogen (ref. 2) is 


Nu b = 0.023 Re®* 8 Pr®* 4 ^ 


\-C, 


( 5 ) 


where 


C = 0. 57 - 

x 

D 


Because of the great interest in hydrogen as a propellant for both nuclear and ad- 
vanced chemical rockets, a large amount of heat-transfer data has been reported. 
These data cover a wide range of conditions, including surface to fluid-bulk tempera- 

ture ratios to 23 and heat flux to 46 megawatts per square meter (28 Btu/sec-in. ). 
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Even though equations (3) to (5) predict friction and heat-transfer coefficients with 
much greater accuracy than has been possible heretofore, some interesting questions 
have been prompted by combining the friction and heat-transfer studies. One concerns 

the possibility of inserting the most recent exponent of T g /T^ from the heat-transfer 
correlation equation (5) into the prior friction correlation equation (4) to make it appli- 
cable to smaller values of x/D. Another point of importance is that the friction equa- 
tion (4) appears to apply to both heat addition (T g /T^ >1.0) and heat extraction 
(V T b < 1.0), but the heat-transfer equation (5) has been tested only for T g /T^ >1.0. 
The question of the applicability of equation (5) to heat extraction (T g /T^ < 1. 0) is cer- 
tainly of interest. This investigation explored and answered these two questions. 


SYMBOLS 


Cj constant used in eq. (6) 

C 2 exponent of T g /T b 

c p specific heat at constant pressure 

D inside diameter of test section 

f/2 half bulk friction coefficient 

fj/2 half film friction coefficient 

G mass flow rate per unit cross-sectional area 

h local heat-transfer coefficient 

k thermal conductivity of gas 

L length of test section 

Nu Nusselt number, hD/k 

Re Reynolds number, VD / v 

Pr Prandtl number, Cpp/k 

T temperature 

V velocity 

X parameter used in reference temperature equation 

x distance from entrance of test section 

i absolute viscosity of gas 

v kinematic viscosity of gas, p/p 

p density of gas 
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Subscripts: 

b fluid bulk (denotes evaluation at fluid bulk temperature, T^) 

cal calculated 

ex experimental 

f film (denotes evaluation at film temperature, = (T g + T^)/2) 

s surface (denotes evaluation at surface temperature, T ) 


CORRELATION EQUATIONS 
Heat -Transfer Coefficients 


Several investigators have used the same basic equation 


Nu 


b 


= c i Re b 8 



( 6 ) 


(with Cj varying from 0. 021 to 0. 024 and C 2 varying with x/D as shown in fig. 2) to 



Figure 2, - Variation of exponent of ratio of surface to fluid- 
bulk temperatures with axial distance from entrance of test 
section. 
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correlate heat-transfer coefficients for helium and air (ref. 3), hydrogen and helium 
(ref. 4), and nitrogen (ref. 5) for T g /T^ to 8. 0. The availability of hydrogen heat- 
transfer data taken over a wide range of conditions led to a thorough study of both the 
data and the many proposed correlation equations (ref. 2). The study recommended the 
use of the correlation equation 


Nu b = 0.023 Re b ,8 Pr b ‘ 4 



( 5 ) 


where 


C 


2 


0. 57 


1. 59 
x 
D 


for all single-phase flow outside the near-critical region and for 2 < x/D < 252. Equa- 
tion (5) was used to calculate more then 3600 heat-transfer coefficients for hydrogen, 
and 87 percent deviated less than ±25 percent from the measured coefficients. In addi- 
tion to the hydrogen data, 88 nitrogen gas and 359 helium gas heat-transfer coefficients 
were predicted with equation (5), and 98 percent deviated less than ±25 percent from the 
measured coefficients. All the data used in reference 2 involved heat addition to the gas 
(T s /T b > 1), and equation (5) was not tested on data with heat extracted from the gas 
Cr s /Tb < 1). 

In the present investigation, heat- transfer coefficients for heat extraction from air and 
carbon dioxide (ref. 6) and air (ref. 7)were used to test equation (5)for its applicability in 
the range 0.25 < T g /T b < 1.0. Predicted heat-transfer coefficients using equation (5) 
are as much as twice the measured coefficients. These results suggest that the term 

^Tg/T^'j 2 doeg not a ppiy to T g /T b <1.0. To determine the effect of T g /T b <1.0 

s 

Nu b = 0. 023 Re 8 ' 8 Pr°* 4 (7) 


on the heat-transfer coefficients, the term [T 




/T. 


b 


was omitted from equation (5) to 


Heat-transfer coefficients predicted using equation (7) are in good agreement with the 
measured data, 81 percent of the calculated coefficients deviating less than ±15 percent 
from the measured values. No effect of T g /T b appears to be present. 
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Reference 6 also reported heat-transfer coefficients for heat addition to air and 
carbon dioxide. Where equation (5) was used to predict these heat-transfer coefficients, 
93 percent of the calculated coefficients deviated less than ±25 percent from the meas- 
ured values. 


Friction Coefficients 


The study of friction coefficients (ref. 1) which recommended the following equation 
for modified surface Reynolds numbers of 3000 or greater 


f 

2 


0. 0007 


0.06251 


Re 


0.32, 



(4) 


used both local and average friction data for both heat addition and extraction (refs. 4 y 
6, and 8 to 12). Equation (4) correlated local friction coefficients for x/D from 
16 to 113. Because equation (5) correlates heat-transfer coefficients for an x/D as 
low as 2 , an obvious step is to use the exponent of T g /T b from the heat-transfer equa- 
tion (5) in the friction equation (4), which gives 


f 

2 


0.0007 



( 8 ) 


where 


D 

Local friction coefficients for heat extraction from air and carbon dioxide (ref. 6) were 
compared with coefficients predicted using equation (8), and the agreement was very 
good for 5 < x/D < 57 (the limit of the experimental data with the exception of the first 
data point which was at an x/D = 1. 57). Equation (8) predicts local friction coefficients 
which are in good agreement with measured values for both heat extraction 
(0. 3 < T g /T b < 1. 0) and heat addition (1. 0 < T g /T b < 7. 35) for x/D from 5 to 113. 
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DISCUSSION OF RESULTS 


Heat -Transfer Coefficients 

The relation shown in equation ( 7 ) was used to calculate heat-transfer coefficients 
for heat extraction from air and carbon dioxide (ref. 6) and air (ref. 7). The coeffi- 
cients for air and carbon dioxide (ref. 6) were calculated using equation (7), with 80 per- 
cent of them deviating less than ±15 percent from the measured coefficients for x/D 
from 5 to 52. Air data (ref. 7) measured only in the entrance section of a tube at x/D 
of 1. 5, 4, 7, and 10 were also calculated using equation (7), with 86 percent of the cal- 
culated values deviating less than ±15 percent from the measured values for x/D of 4, 

7, and 10. The calculated coefficients at x/D of 1. 5 ran 15 to 30 percent lower than 
measured values. 

The additional data for heat added to air and carbon dioxide (ref. 6) were local 
heat-transfer coefficients for x/D from 0.22 to 59 and Tg/T^ to 2. 8. Using equa- 
tion (5), 98 percent of the calculated heat-transfer coefficients for air deviated less than 
±15 percent for x/D from about 2 to 59. The calculated values for carbon dioxide for 
the same x/D did not agree as well as the air data, with 88 percent deviating less than 
±25 percent. Table I summarizes the results of both the present investigation and ref- 
erence 2. 


TABLE I. - PERCENT OF PREDICTED LOCAL HEAT -TRANSFER 
COEFFICIENTS THAT DEVIATE LESS THAN ±25 PERCENT 
FROM EXPERIMENTAL COEFFICIENTS 


FOR VARIOUS GASES 


Gas 

Number of data points 

Percent of h ca ^ that deviate less 

than ±25 percent from h 

1 ex 

Heat added to gas, 1. 0 < T g /T^ < 23; calculations used eq. (5) 

Hydrogen 

3674 

87 

Helium 

359 

98 

Nitrogen 

88 

97 

Air 

370 

97 

Carbon dioxide 

370 

88 

Heat extracted from gas, 0. 25 < T g /T 

k < 1.0; calculations used eq. (7) 

Air 

641 

86 

Carbon dioxide 

550 

94 
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Friction Coefficients 


The relation shown in equation (8) was used to calculate 1100 friction coefficients 

for heat extraction from air and carbon dioxide (ref. 6) for 0. 3 < T /T^ <1*0 and 
5 < x/D < 52. Of the 1100 friction coefficients calculated, 94 percent deviated less than 

±15 percent and 88 percent deviated less than ±10 percent from the measured values. 

Figure 3 shows the product of the friction coefficient and the temperature ratio cor- 
rection factor as a function of the modified surface Reynolds number. A total of 1523 


Fluid 


O Air 
Air 
Q Air 
# Air 


Type of friction coefficient 

Local and average 
Local 
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Reference 
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Unpublished data of Wetland (of Lewis) 
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□ Helium Average 
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Figure 3. - Correlation of local and average friction coefficients for modified surface Reynolds 
numbers of 300Q and greater. Density in friction coefficients evaluated at bulk temperature; 
viscosity in Reynolds number evaluated at surface temperature; 1523 data points in turbulent 
region. In ordinate D/x applies to local coefficients, D/L should be used for average coefficients. 


friction coefficients for both the heating and cooling of gases were calculated using 
equation (9); 95 percent of the data are within ±15 percent and 89 percent fall within 
±10 percent of the correlation line. The data in figure 3 indicate a lack of a transition 
region, except for the few points for air with a T g /T b <1.3 that fall below the corre- 
lation line. Thus, it appears that the large variations in the fluid properties tend to 
make the flow turbulent at modified surface Reynolds numbers of 3000 or more. 
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CONCLUSIONS 


Single-phase heat-transfer and friction coefficients for hydrogen, helium, nitrogen, 

air, and carbon dioxide with large variations in the physical properties have been 
studied. The coefficients were measured over a wide range of conditions outside the 

near-critical region for subsonic flow through smooth straight tubes. For the condi- 
tions stated, the following conclusions may be made: 

1. For heat extraction from air and carbon dioxide with 0. 25 < T /T. <1.0 and 
4 < x/D < 52, there appears to be no discernable effect of either x/D or T g /T b on 
the heat-transfer coefficient and the equation 

Nu b = 0. 023 Re®' 8 Pr b * 4 

where Nu fe is the bulk Nusselt number, Re b is the bulk Reynolds number VD /u^, and 
Pr b is the bulk Prandtl number, is recommended to predict heat-transfer coefficients. 
Using this equation, 80 percent of the calculated heat-transfer coefficients deviated less 
than ±15 percent and 90 percent deviated less than ±25 percent. 

2. For heat addition to hydrogen, helium, nitrogen, air, and carbon dioxide with 
1. 0 < T g /T b < 23 and 2 < x/D 252, there is a strong effect of both x/D and T g /T b> 

A total of 4900 heat-transfer coefficients were calculated with the relation 

Nu b = 0.023 Re»- 8 Pr»' 4 ^j 

where 

C 9 = 0. 57 - 

Z X_ 

D 

and Nu b , Re b , and Pr b are the bulk Nusselt, Reynolds, and Prandtl numbers, re- 
spectively,- T g is the surface temperature, T b is the fluid-bulk temperature, x is the 
axial distance from the entrance, and D is the inside diameter of the tube. Ninety per- 
cent of the calculated coefficients deviated less than ±25 percent from the measured 

values. 

3. Friction coefficients for heat addition to hydrogen, helium, nitrogen, air, and 
carbon dioxide with T g /T b to 7. 35 and 5 < x/D < 200 and heat extraction from air 
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and carbon dioxide with T g /T b as low as 0. 30 and 5 < x/D < 57 all show an effect of 

both x/D and T /T, and have been accurately predicted by the relation 
S D 



D 


where Re is the modified surface Reynolds number VD/n , T_ is the surface tern- 
perature, is the fluid-bulk temperature, x is the axial distance from the entrance, 
and D is the inside diameter of the tube. Of the 1100 friction coefficients calculated, 
88 percent deviated less than ±10 percent and 94 percent deviated less than ±15 percent 
from the measured values. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, December 30, 1970, 
122-29. 
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